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Abstract 

We propose a natural expansion of the atomic field operator in 
studying elementary excitations in trapped Bose-Einstein Condensa- 
tion (BEC) system near T=OK. Based on this expansion, a system of 
coupled equations for elementary excitations, which is equivalent to 
the standard linearized CP equation, is given to describe the collective 
excitation of BEC in a natual way. Applications of the new formalism 
to the homogeneous case emphasize on the zero mode and its relevant 
ground state of BEC. 
PACS number(s): 03.75.Fi, 05.30.Jp 



1 Introduction 

Since the realization of Bose-Einstein Condensation (BEC) in trapped alkali 
atomic vapors in 1995 much attentions have been attracted from both 
viewpoints of theories and experiments Q. Among all the recent develop- 
ments, the study of collective excitations occupies an important position in 
studying the properties of trapped BEC. When the temperature approaches 
the critical temperature, the theoretical description of collective excitation 
is still proceeding at present |^]. For the case that the temperature becomes 
near zero, the linearized GP equations ||7|, |8[ can give correct numerical pre- 
dictions in agreement with the experiments jl^ . However, the linearized GP 
equations are difficult to be solved analytically, and they include confused 
zero mode problemQ. In order to give some new light on this problem, we 
develop an equivalent formalism to the linearized GP equations based on a 
natural expansion of the atomic field operator. 



Our study follows discussions in ref.[P by Lewenstein and You, and the 
main difference is that we adopt specific eigenfunctions to expand the atomic 
field. It is noticed that the advantages of our expansion lie in the follow- 
ing facts: (1) The problem in solving the elementary excitation becomes 
the problem in diagonalizing the standard quadratic Hamiltonian of Boson 
operators, which has been completely solved in principle before^. (2) Our 
formalism adopts a simpler form than those formulated by other expansion 
of the atomic field operator. (3) Since only the several lowest excitations 
dominate the behavior of trapped BEG system near zero temperature, the 
limited-level approximation (which is explained in details in Sec. 2) is ra- 
tional in most cases. 

This paper is organized as the following: A general formalism for ele- 
mentary excitations is present in Sec. 2, and then in Sec. 3 it is applied to 
the homogeneous case, which is interest since it can be used to be a good 
example to demonstrate the appearance and the physical meaning of the 
zero mode clearly. In Sec. 4, a conclusion will be made and some possible 
application of our special expansion will be discussed. 



2 General Formalism 

The Hamiltonian of trapped BEG system is 

r - h'^v'^ ^ 1 . . . . 

H = y d3f[^H^)(-^^ + ^tr(r)-Mo)V'(0 + 25^H^)V'H^)V'(^=')^(^)], (1) 

where /io is the chemical potential, g = , Ug the s-wave scattering 

length of the interatomic potential, M the atomic mass, Vtr{r) the external 
potential, and ip{r) ('(/'^(r)) the atomic annihilation (creation) field operator. 

Near T = OK, since almost all the atoms occupy the same one-particle 
quantum state ^o{'r), it is convenient to separate out the special mode for 
the ground state of BEG from the atomic field operator 

V'(r) = yiVo$o(r) + 5V^(f). (2) 

where 5tl^{f) represents the quantum fluctuation of tp^f) relative to\/ No^o^f). 
Notice that, at present J d^r (Sip'^ {r)5ip{r)\ <^ Nq, i.e. the thermal 

\ / ensemble 

component atoms is negligible. Hence the Hamiltonian can be properly ex- 
panded in 6Tp{'r) series while substituting formula(^) into formula(|l|) . (Only 
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the lower order of 6ilj{r) is important, so up to two orders of 5ij){r) is main- 
tained in our consideration.) From the above Hamiltonian, ^^{r) can be de- 
termined in zero order of 5tp{r) by making use of variation methods, which 
satisfy the time- independent GP equation [P, |^ 

^" W + ^^^^''^ + gNQ%{r)^^{m^{r) = /Uo<&o(r). (3) 

Combining eqs. (IH-^), and ignoring the terms of the third and fourth 
order of 6ijj{f) in eq.(|T|), the Hamiltonian can be simplified as. 



+^gNo J dV#t(f)5V;t(f)$o(r)^o(^. (4) 

In order to further simplify the Hamiltonian, we select a set of special 
complete wave functions {<^„(r)} to expand the atomic field operator, 

'^i^) = Yl ^n^nir) + Ao<!>o{f), (5) 
where{<I>„(r')} satisfy the following equations 

^" W + ^''^"^ + <7iVo^S(r1$o(r))1>n(r) = /i„$„(r-), (6) 

Un is the annihilation boson operator of the single state ^nif') {n 7^ 0), Aq 
the annihilation boson operator of the single state ^^{r). Due to the fact 
that (--^^ + Vtrir) + gNQ(^l{r)(tQ{r)) is Hermitian, its eigenstates <I>„(f) 
(n = 0, 1, • • •) are orthorgonal and complete 



dV$:;,(f)$„(r) = 5^mn. (7) 

In most cases, only the lower part of elementary excitations determines the 
properties of the system near zero temperature, so it is a good approxima- 
tion to hold only the lower /-level of eq.@, i.e. n = {0, 1, •••,/ — 1}. In 
addition, when / 00, the present treatment becomes rigorous in principle. 
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Substituting eq. (^) into eq. (^), we obtain 

/-I 

5V^(r) = ^aA(f), (8) 

n=0 

where 

ao = io - v^. (9) 
In terms of eqs. the Hamiltonian (^) can be reexpressed as 

m.n 



2 



where 



5mn = ffiVo y dV-cI>o(f)cI>o(f)cI>;;(f)c^;(r), (12) 

dmn = gNo j <fr^U^%{^^o{f^^n{r). (13) 

The scheme of diagonahzing the Hamiltonian as the form ([To|) has been 
extensively studied by Blaizot and Ripka[Q] . Here, their results will be briefly 
reviewed in the following for the use in our present discussion. 

First, one write the Hamiltonian in a compact form 

H ~ -^^^0 + ]^Q^Ma - hrA (14) 
by introducing the vector operator 

= ( a"!" a ) = ^ Oq a\ ■ ■ ■ ao " " " ^ , (15) 

the 2/ X 2/ coefficients matrix 

^={b* a* )'^=^^^ = ^• (16) 

In ref.0], it is assumed that the matrix M is semi-definite positive. However, 
in our discussion, this condition on the matrix M is not necessary since it is 
only a condition for the stability of the mode of the ground state of BEC M . 
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Second, a unitary canonical transformation is carried out to diagonalize 
the Hamiltonian, 

p = Ta (17) 

where the operator vector 

f3^ = {b^ b) = [bl hi ■■■ b\_^ bo h ■■■ (18) 
the transformation matrix and its inverse matrix 

the matrix satisfy the canoical condition 

Tt/T^? = 1, (20) 

and the matrix 

If the elements of the / x / matrixes X and Y satisfy the following 
conditions 0, 14] 

i]MV'' = WnV'\ (22) 

where 

f l,xr = x„„y," = y„„(z = i,2,---,/) 



especially, if the above equation have one special solution {wq = 0,V^ = P}, 
i.e. 

r]MP = (23) 
the Hamiltonian can be written as 

H :^ -^No + uJib^ + g + i ^ ^„ _ ItrA (24) 

71=1 n 

where 

p = a^r]P, (25) 

and ^ is a positive constant, which can be determined by the following 
conditions 

r/MQ = -i-, (26) 

Q^VP = i, (27) 
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where the vector Q is orthorgonal to ah the eigenvectors of the matrix rjM. 

In Hamiltonian (|2^), the second term is the Hamiltonian of a system of 
independent oscillators, which represent elementary excitations of the sys- 
tem; However, the third term has the form of a free kinetic energy, which 
is connected with a collective motion in Fock space arising from a broken 
U{1) symmetry in the procedure of the mean field approximation Usu- 
ally, the third term is termed with spurious state^ or zero mode[^ for the 
corresponding eigenvalue and the norm of the vector P both are zero. The 
physical meaning of this term will be discussed in details in Sec. 3. 

In fact, we can determine the vector P through observation. Notice that 
in principle, we can solve the the eigenvalues {cOn} and the corresponding 
eigenvectors defined by {X^,y^} of eqs. (p^, which consist of 2 x / x / 
homogeneous linear equations. Obviously, the eigenvalue = and the 
corresponding vector P denoted by = 6mO,Y^ = —6mo] is a specific 

solution of the above equations. Hence 

p = ao + 4; (28) 

which is in agreement with that in ref.[Q. 

From eqs. ( p4| , p^ ), a direct conclusion is that the approximate vacuum 
state \Vac) = Hi ® 1^^^)^ of BEG satisfies the following conditions 

k\Vac)^ = 0, (29) 
p\Vac)^ = 0. (30) 

To sum up, in this section, we give a new formalism for elementary 
excitations in trapped BEG, which is equivalent to the standard linearized 
GP equation. This equivalence can be easily verified when the complex wave 

(uiv) \ 
,^ in traditional method are expanded with the specific 
v{r) ) 

complete wave functions defined by eqs.(P). In this sense, our formalism is 
a specific representation of the traditional method. 

3 Homogeneous case Vtrif) =0 

In this section, the general formalism obtained in the above section will be 
demonstrated in the homogeneous case. 

In this case, the ground wave function satisfy 

^" W + 5iVo^S(^"')^o(f))1>o(r) = ^o$o(r). (31) 
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Therefore, the ground wave function and the chemical potential are given 
by 

$o(r) = 47'M = 5^. (32) 



(33) 



ly " - V 

Eq. (0) which give the complete wave functions now becomes 



(- 



' 2M 



By solving the above equations, the eigen wave functions and the corre- 
sponding eigen values are given by 



1 



Ak-r 



2 1,2 



2M 



(34) 



Formulas (P,IT3D become 



^kk' 



^kk' 



gNo I dV<I>o(f)cI>o(f)cI>t(f)cI,t,(f)=5^%_fc,, 



gNo I dVc^.t(f)$S(f)$o(f)$,^,(f)=5^%,. 



(35) 



Due to the fact the mode denoted by k is only coupled to the mode denoted 
by —k implied by the above equation, eqs.(|2^) can be simplified as 



h^k^ 



No 

2M ^ V 

^2^.2 



2M 



V 



(36) 



When A;' 7^ 0, the eigenvalue can be calculated by requiring that the 
above equations have nontrivial solution, 



h^k"^ 



+ 9 
2M ^ V 



0\2 



(37) 



The corresponding annihilation operators of the elementary excitation are 



^(^^^^^ + l)a„ 



-( - l)aK 
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or 



2M ~^ y V 



+ l)a 



1 Mfc^ _L 

^1 2M ~^ y V 
2^' 



IjJt 



Clearly, it comes back to the familiar form|^], which supports that our for- 
malism is equivalent to the traditional one. 

Since the operators do and oj are only coupled each other, we can limit 
ourself in the subspace of the wave vector A;' = 0. The matrix r/M in this 
subspace is 



^M 



y V y V 
y V y V 



Obviously, the eigen vector P of the zero mode and the corresponding mo- 
mentum operator p are obtained as 



P 



1 



p = ao + a^. = {Ao + A^) - 2 ViVo- 



According to eqs. (26,27), the other independent vector Q and the constant 
/i are obtained as 



Q = — 

^ 2 



In sum, the Hamiltonian can be written as 

2// 



2V ^ *:fc^ 2ii 2^ 



2 Z d k 

fc^O 



No 



Now, the approximate vacuum state l^ac) = Hfc ® "-"^ BEC can 

be obtained analytically by solving eqs.( p9| , ^ ). When the wave vector k ^ 
0, the vacuum state \Vac)^ is given by solving the equations 6^|yac)g = 
hAVac)r = Q, 



2M "T y y 



2M ^ y V ^^k, 



\n)r (g) |n)_ 



where Ar is the constant of normalization, the vector |n) 



■k ' 



is the 



eigenstate of the number operator dgd^ (d' When the wave vector 

k = 0, p \ Vac)Q = 0. Therefore, 
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\Vac), = -^Y.\^)^l dxe^^^^x\nr, (39) 
where the state vector |n)Q is the eigen state of the number operator Aq^Qj 

and Hn{x) is Hermit polynomial. 

In this section, we explicitly solve the elementary excitations of BEC 
in the homogeneous case. In this case, it is easy to see that the term of 
the zero mode in Hamiltonian ( |38| ) originate from the quantum fluctuation 
of the mode denoted by macroscopic wave function ^o{f). In fact, if we 
adopt the usual Bogoliubov approximation, i.e. V^o^ thus the 

momentum operator p = 0. However, we have no specific reasons to ignore 
this quantum fluctuation while maintaining those of the other modes. In 
our treatment, the conservation of the particle number is destroyed, which 
is easily seen from the approximate vacuum state l^ac). The kinetic term 
appearing in the Hamiltonian is originated from this symmetry breaking, 
which represents a collective motion, not an intrinsic elementary excitation 
of the system Q. 



4 Conclusion 

In this paper, based on a natural choice of the complete wave functions, 
we expand the atomic field operator and obtain a new formalism for the 
excitations of trapped BEC system near zero temperature. We argue that 
our formalism is equivalent to the standard linearized GP equation. In 
terms of this formalism, we illustrate the relation between the zero mode 
and the other excited modes. Essentially, the zero mode originates from the 
quantum fluctuations of the mode denoted by the condensate wave function. 
When applicating the formalism to the homogeneous case, the formalism 
comes back to the usual Bogoliubov excitation spectrum, which identifies 
our theory.. Especially, in this case, the physical meaning of zero mode 
become obvious and the ground state of BEC can be calculate explicitly up 
to second order of the quantum fiuctuations. 
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